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Abstract. - Heat fluctuations of a harmonic oscillator in contact with a thermostat and 
driven out of equilibrium by an external deterministic force are studied experimentally and 
theoretically within the context of Fluctuation Theorems. We consider the case of a periodic 
forcing of the oscillator, and we calculate the analytic probability density function of heat 
fluctuations. The limit of large time is discussed and we show that heat fluctuations satisfy the 
conventional fluctuation theorem, even if a different fluctuation relation exists for this quantity. 
Experimental results are also given for a transient state. 



Out-of-equilibrium fluctuations have recently received a lot of attention, especially in the 
context of nanosystems and biomolecules where fluctuations are large. In the last decade, 
Fluctuation Theorems (FTs) appeared in noncquilibrium physics. They quantify the asym- 
metry of fluctuations of entropy production for a large class of systems, possibly far from 
equilibrium. These theorems were first demonstrated in deterministic dynamics of many de- 
grees of freedom [1,2] and later extended to stochastic systems [3-6]. The FT states that the 
probability P(a = a) of observing an entropy production rate a, measured over a time r, with 
a value a, satisfies 

-p^ — > exp(ar) for large r and any value a (1) 

There are two classes of FTs. The Stationary State Fluctuation Theorem (SSFT) considers 
a noncquilibrium stationary state. SSFT holds in the limit of large time r. The Transient 
Fluctuation Theorem (TFT) describes transient noncquilibrium states where r measures the 
time interval since the system left its equilibrium state. Contrary to SSFT, the TFT holds 
for all times, i.e., equation JQl is an equality even for arbitrarily small values of r. 

In this letter, we consider the heat Q T dissipated on a time interval r in a thermostated 
system at temperature T rather than the entropy production. We study the fluctuations 
of the heat dissipated by a harmonic oscillator in contact with the thermostat and driven 
out of equilibrium by an external force. Experimental tests of FTs are rare. Some of them 
are performed in dynamical systems [8], in which the interpretation of the results is very 
difficult. In stochastic systems, a laboratory experiment was carried out by Wang et al using 
a Brownian particle in a moving optical trap [9]. Work fluctuations were shown to obey the 
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predictions of ref. [5] . Work and heat fluctuations were also studied in an electrical circuit by 
Gamier and Ciliberto [10] ; the theoretical predictions for both heat and work fluctuations 
were satisfied [5,6]. These two systems are described by a first order Langevin equation. 
Systems described by a second order Langevin equation have been studied [3], and tested 
experimentally for work fluctuations W T [14]; as far as we know, no analytical results for the 
probability density functions (PDFs) of heat fluctuations Q T are available in this case. In 
ref. [5] , van Zon and Cohen have calculated the Fourier transform of the PDF of Q T for a first 
order Langevin dynamics but no exact expression of the PDF itself is known. 

In the following, we consider first a stationary state obtained by driving the system period- 
ically in time. We calculate exactly the probability density function of Q T . We then compare 
our finding with new experimental results, and show that SSFT holds, and we discuss the 
large time limit. Finally experimental results of a TFT for the heat are reported. 

To study FT, we measure the out-of-equilibrium fluctuations of a harmonic oscillator 
whose damping is due to the viscosity of a surrounding fluid, acting as a thermal bath. The 
oscillator is a torsion pendulum composed of a brass wire and a glass mirror glued in the 
middle of this wire. The elastic torsional stiffness of the wire is C = 4.7 x 10~ 4 Nmrad -1 . 
It is enclosed in a cell filled by a water-glycerol solution at 60% concentration. The system is 
a harmonic oscillator with resonant frequency f a = \J C / I e g / (2-k) = 217 Hz and a relaxation 
time r a — 2I c ff/v — 9.5 ms. I e s is the total moment of inertia of the displaced masses and v 
is the oscillator viscous damping. The angular displacement of the pendulum 9 is measured 
by a differential interferometer [12, 13]. The measurement noise is two orders of magnitude 
smaller than thermal fluctuations of the pendulum. 9{t) is acquired with a resolution of 24 
bits at a sampling rate of 8192 Hz, which is about 40 times f . The calibration accuracy of 
the apparatus, tested using the Fluctuation Dissipation Theorem, is better than 3% (see [13]). 
We drive the system out-of-cquilibrium by forcing the system with an external torque M by 
means of a small electric current J flowing in a coil glued behind the mirror. The coil is inside 
a static magnetic field, hence M oc J. More details on the set-up can be found in ref. [12,13]. 
The system is very well described by a second order Langevin equation: 



with r/ the noise, delta-correlated in time, /3 = (fc^T)^ 1 with ks Boltzmann's constant and 
T the temperature of the system. At equilibrium (M = pN.m), the PDF of the thermal 
fluctuations SO is Gaussian with variance ksT/C. We apply two kinds of forcing. First, we 
study a periodic forcing : M(t) — M Q sm(ojdt) with M a — 0.78 pN.m and Wd/27r = 64 Hz. We 
then analyze a linear forcing M(t) = M t/r r with M Q — 6.2 pN.m and r r = 0.01 s = 1.07 r a . 

The change in internal energy AU T of the oscillator over a time t, starting at a time ij, is 
written as: 



which is the first law of thermodynamics. W T is the work done on the system over a time r: 



d 2 d9 



+ C6 = M + ^2k B Tur], 



(2) 



AU T = U(U + r) - U{U) = Q T + W T 



(3) 




(4) 



and Q T is the heat given to the system. Equivalently, (— Q T ) is the heat dissipated by the 
system. 

For a harmonic oscillator described by a second order Langevin equation, the internal 
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energy has two contributions : the kinetic and potential energies: 



U(t) 



1 



k B T 



1 



d0(t) 



dt 



1 



C8(t) 2 



(5) 



Multiplying eq. (J2J) by 41 and integrating between ti and ti + r, we obtain exactly the first 
law of thermodynamics eq. and have the following expression for the heat: 



Q T = AU T -W r = - 



k R T 



ti+T 



dt' 



k H T 



U+T V(t')^(t')dt' . (6) 



it 



The first term corresponds to the viscous dissipation and is always positive, whereas the second 
term can be interpreted as the work of the thermal noise which have a fluctuating sign. 

We rescale the work W T (the heat Q T ) by the average work (W T ) (the average heat (Q T )) 
and define: w T — (q T — Jq~))- Averages are obtained experimentally as time-averages, 

and they are proportional to r on the stationary state under consideration. 

We consider first the periodic forcing. In this case, we choose r as a multiple of the period 
of the driving (r = 2mr/ujd with n integer). The starting phase Uujd is averaged over all 
possible ti to increase statistics. The PDFs of w T , AU T and q T are plotted in Fig. for 
different values of n. The average of AU T is clearly zero because the time r is a multiple of 
the period of the forcing. The PDFs of the work (fig. are Gaussian for any n whereas the 
PDFs of heat fluctuations q T have exponential tails (fig. ^). These exponential PDFs can be 
understood noticing that, from eq. 10, — Q T = W T — AU T and that AU T has an exponential 
PDF independent of n (fig. 2Jd). Therefore, on a first approximation, the PDF of q T is the 
convolution between an exponential and a Gaussian. 

To quantify the symmetry of the PDF around the origin, we define the function S as: 



S{e T ) 



1 



(Er) 



hi 



P(e T ) 



P(-er) 



(7) 



where e T stands for either w T or q T and E T stands for either W T or Q T . The question we ask 
is whether: 

lim S(e T ) = e T (8) 

r — >oo 

as required by SSFT. S(q T ) is plotted in Fig.^i for different values of n ; three regions appear: 

(I) For large fluctuations q T , S(q T ) equals 2. When r tends to infinity, this region spans 
from q T = 3 to infinity. 

(II) For small fluctuations q T , S(q T ) is a linear function of q T . We then define S g (n) as the 
slope of the function S(q T ), i.e. S(q T ) = T, q (n)q T . This slope is plotted in Fig. ^ where we 
see that it tends to 1 when r is increased. So, SSFT holds in this region II which spans from 
q T = up to q T = 1 for large t. 

(III) A smooth connection between the two behaviors. 

The PDF of the work being Gaussian, the functions S(w T ) are proportional to w T for any r, 
i.e. S(w T ) = Yj w (n)w T (ref. [14]). T, w (n) is plotted in Fig.^3 and we observe that it matches 
experimentally S 9 (n), for all values of n. So the finite time corrections to the FT for the heat 
are the same than the ones of FT for work [14] : T, w (n) = ~S q (n) = l + K/n + 1/nO (e -r "/ T ") , 
where if is a constant. 

We now give an analytical expression of the PDF of heat Q T . In order to do this, we 
write 9 = 9 + 89 where 9 is the mean response of the system to the external work and 59 the 
thermal fluctuations. We suppose that fluctuations at M(t) ^ are those of equilibrium, i.e. 
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w T = W T /<¥ T > 
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q T =Q T /<Q T > 




Fig. 1 - Sinusoidal forcing, a) PDFs of w T . b) PDFs of AU T . c) PDFs of q T . d) Functions S(q T ). 
In all plots, the integration time r is a multiple of the period of forcing, r = 2nn/uid, with n = 7 
(o), n = 15 (□), n — 25 (o) and n = 50 (x). Continuous lines are theoretical predictions with no 
adjustable parameter, e) The slope £<j(n) of S(q T ) for q T < 1, plotted as a function of n (o). The 
slope T, w (n) of S(w T ) plotted as a function of n(D). Continuous line is theoretical prediction. 

that the external driving does not perturb the equilibrium PDF. This hypothesis is supported 
experimentally as shown in ref. [13] and detailed in [14]. Times r under consideration are taken 
larger than r Q = 9, 5 ms, so that exponential corrections to the autocorrelation functions, 
which are scaling like exp(— r/r Q ), can be neglected. Experimentally, r/r a = 1.64 n, so this 
is correct as soon as n is larger than 3 or 4. Within this assumption, 6(ti + r) and 6(ti) 
are independent, and so are 4f(ij + t) and 4| (ti). Eq. J2| is second order in time, so 6(t) 
and 4j(t) are independent at any given time t. Just like in our experiment, we choose the 
integration time r to be a multiple of the period of the forcing, so (AU T ) = and therefore 
(W T ) = —(Q T )- Within this framework, we find that the PDF of Af7 r is exponential: 



P(A£/ T ) = 7.exp(-|A<7 r |) 



(9) 



It is independent of r because AU T depends only on 9 and 4j at times U and U + t which 
are uncorrelated. This expression is in perfect agreement with the experimental PDFs for 
all times (see Fig. \3p). Experimentally, work fluctuations have a Gaussian distribution so it 
is fully characterized by its mean (W T ) and its variance a^- In a Gaussian case, the FTs 
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take a simple form S(w T ) = 2 2 w T = 'E w (n)w T . We have already computed the analytic 

expression of H w {n) [14] and this gives a relation between the variance and the mean value of 
W T : o 2 w = 2(W T ) + O(i) and so a 2 w = 2\(Q T )\ + 0(1). 

To obtain the PDF P(Q T ) of the heat, we define its Fourier transform (characteristic 
function) as 

/>OC 

P T {s) = / dq T e is ^P(q T ) (10) 



which can be computed exactly [5]. We then write P(q T ) using eq. JBJ as: 

P{q T )= [ fd069P(AU T -Q T ,O(U+T),e{U)AU+r),e(ti)) (H) 



where P is the joint distribution of the work W T , 6 and 4f at the beginning and at the end 
of the time interval r. This distribution is expected to be Gaussian because W T is linear in 9 
and additionally 9, 9 and W T are Gaussian. Some algebra then yields: 



1 , a 



2 



PA*) = ^cxp \j(Q r )s - -fs 2 j (12) 

The characteristic function of heat fluctuations is therefore the product of the characteris- 
tic function of an exponential distribution (■jxpO with the one of a Gaussian distribution 

(exp (i(Q T }s — ^-s 2 ^j). Thus the PDF of heat fluctuations is nothing but the convolution of 

a Gaussian and an exponential PDF, just as if W T and AU T were independent. The inverse 
Fourier transform can be computed exactly: 



3 Q.-(Q.) erfc ( 9jL 




(Qr-(Qr)) eTic Z*l Z WW I » jy f (13) 



where erfc stands for the complementary Erf function. In Fig. ^j, we have plotted the ana- 
lytical PDF from eq. (|13|l together with the experimental ones, using values of and (Q T ) 
from the experiment. The agreement is perfect for all values of n and with no adjustable 
parameters. Using eq. (fl^jl . we isolate three different regions for S(q T ): 

(I) if Q T > a\ r + \ {Q T )\ = 3|(Q T )| + 0(1), then S(q T ) = 2 + O(i). This case corresponds 
to large fluctuations and the PDF can be pictured as exponential with a non-vanishing mean. 

(II) if Q T <a 2 w ~ \(Q T )\ = \(Q T )\ + 0(1), then S(q T ) = S(n)g T + 0(±) with E,(n) = 
T ^ = H w (n). In this case, heat fluctuations are small and behave like work fluctuations. 

The slope S(r) is the same as the one found in the case of work fluctuations. The exact 
correction to the asymptotic value is plotted in Fig. ^3 and again it matches perfectly the 
experimental behavior. 

(III) for afy — \ (Q T ) \ < Q T < Oyv + KQt)|) there is an intermediate region connecting cases 
(I) and (II) by a second order polynomial S(q T ) = 2. - ^-{q T - (1 + sff))) 2 + 

Those three regions offer a perfect description of the three domains observed experimentally 

(Fig.nn). 

Now, we examine the limit of infinite r in which SSFT is supposed to hold but which 
depends on the variable we use : either the heat Q T or the normalized heat q T . First we discuss 
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Q T . The asymptotic behavior of the PDF of Q T (eq. H13|) ~) for large r is Gaussian with variance 
a^y. Thus, the PDF of Q T coincides with the PDF of W T for r strictly infinite. As we have 
already shown, work fluctuations satisfy the conventional SSFT, therefore heat fluctuations 
also satisfy the conventional SSFT (eq. Q). We have found three different regions defined 
by two characteristic values, but in the limit of infinite time r, only region (II) is relevant: 
region (II) is bounded from above by |(Qt)| + 0(1) with the average (Q T ) being linear in 
r. We see that all the behavior of the fluctuations of Q T is captured by region (II) where 
S(q T ) is linear for all Q T and SSFT holds. Second, we turn to the normalized heat q T . As 
the average value of Q T is linear with r, rescaling by (Q T ) is equivalent to a division by r ; 
the mean of q T is then 1. This normalization changes into constants the two characteristic 
values which delimit the three regions : the boundary between (II) and (III) is now l + C(l/r) 
and the boundary between (III) and (I) is 3 + C(l/r). The function S(q T ) is not linear for 
large values of q T but only in region (II) (q T < 1), for small fluctuations. So SSFT is satisfied 
only for small fluctuations but not for all values of q T , and instead of a FT, we rather have 
a fluctuation relation. These two descriptions in terms of Q T or q T are in fact a problem of 
two non-commutative limits. The first description implies that one takes the limit r infinite 
before taking the limit of large Q T . The second description does the opposite. However, the 
probability to have large fluctuations decreases with t and experimentally, for large r, only 
the region (II) where SSFT holds (q T < 1) can be seen. 

Finally, we briefly report results for the Transient Fluctuation Theorem. For this, we 
choose a torque M(t) — M t/r r linear in time, and the system is at equilibrium at tj = 
(M(ti = 0) = pN.m). Unlike with the periodic driving, the average of the variation of 
internal energy (AU T ) is not vanishing. The work done by M(t) is used by the system to 
increase his internal energy but a small amount of energy is lost by viscous dissipation and 
exchange with the thermostat. The PDF of W t /t, AU t /t and Q t /t are plotted in Fig. [3 
for different values of r/r a . Averages (W t /t) and (AU t /t) are linear in r. However, their 
difference (eq. (0) (Q t /t) is constant and of the order of a few ksT/s (Fig. [2:). The shapes 
of the PDF are different from the ones obtained with the periodic forcing. Work fluctuations 
have a Gaussian PDF for any values of r, moreover TFT holds for W T [14]. In contrast 
AU T have a different probability distribution. Fig. 03 shows that the PDF are not symmetric 
around the mean value. Extreme events have again an exponential distribution. For exactly 
the same reasons as in the case of a periodic forcing, the PDFs of Q t /t are not Gaussian and 
have exponential tails for extreme fluctuations : P(Q t /t) — Aexp(— a|Q T /r|). As we can 
see in fig. the PDFs are not symmetric around the mean value, thus there are two pairs 
(a, A), one for positive value of large fluctuations (a+, A + ) and one for negative («-, ^4-)- 
Thus, there is a simple expression of S(q T ) for large fluctuations: 



It can be seen experimentally in fig. S(q T ) is not proportional to q T , therefore TFT is 
not satisfied for finite time. However, for large value of r, the PDF of Q T becomes symmetric 
around the mean value and only the Gaussian behavior remains. Thus, TFT appears to be 
satisfied experimentally in the limit of infinite r. This breaks the expected property of TFT 
to be valid at any time. 

In conclusion we have proposed an SSFT for heat fluctuations of a harmonic oscillator 
driven in a stationary out of equilibrium state by a periodic external force. An exact expres- 
sion of the PDFs of the heat Q T averaged over a time r is given. This PDF is in perfect 
agreement with experimental data. For finite times, we have isolated different behaviors: one 




(14) 
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lQ T (k B T/s) q T =Q T /<Q T > 



Fig. 2 - TFT. All plots are given for two values of r/r Q : 0.3 (□) and 1.04 (x). a) PDFs of \W T . 
Continuous lines are Gaussian fits, b) PDFs of -£JJ r - c) PDFs of -Q T . Continuous lines are 
exponential fits for large Q T . d) Functions S(q T ) for large heat fluctuations q r computed using 
exponential fits. 

for small fluctuations (Gaussian behavior) and the other for extreme fluctuations (exponen- 
tial behavior). SSFT holds for infinite time. We have also studied a TFT for Q T using linear 
forcing and found that FT is satisfied only in the limit of large times. 

We thank G. Gallavotti for useful discussions. This work has been partially supported by 
ANR-05-BLAN-0105-01. 
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